Initial orientation effect and selecting desired events 
in 520AMeV/u U-U collisions 



K.J. Wu,i'2 F. Liu,i'2 and N. Xu^ 

^ Institute of Particle Physics, Hua-Zhong Normal University, Wuhan 430079, China 
^ Key Laboratory of Quark & Lepton Physics (Huzhong Normal University), Ministy of Education, China 
Nuclear Science Division, Lawrence Bekeley National Laboratory, Bekeley, CA 94720, USA 

How to select out those collisions with the desired geometry such as tip-tip and/or body-body in 
experiment is one key point for performing high energy UU collisions. With a relativistic transport 
model, we performed a simulation for deformed UU collision with vast different orientations at 
CSR energy area corresponding to the high net-baryon density region in QCD phase diagram. By 
investigating the centrality and initial collision orientation dependence of the center baryon density, 
we found that the tip-tip like UU collisions with extended high density phase, which is very important 
for studying the nuclear EoS of high baryon density matter and the possible end-point of the phase 
boundary, are those events with small initial orientations (< 20°) for bath projectile and target 
in reaction plane and small impact parameter (< 2.6/m). We pointed out quantificationally two 
observations - multiplicity of forward neutron and nuclear stopping power that both allows us to 
select out those most interesting events (i.e. tip-tip like), which will be very helpful for the future 
experiments at performing UU collisions. 
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I. INTRODUCTION 

A major goal of current and future high energy heavy- 
ion collisions experiments is to probe and study the state 
of new matter under extreme conditions. Lattice QCD 
calculation predicts that a phase transition from the 
hadronic state to the quark-gluon-plasma (QGP) will oc- 
cur under high temperature or density. 

In recent years, the focus on searching for QGP is 
performed at SPS/CERN and RHIC/BNL with ultra- 
relativistic heavy ion collision correspond to high tem- 
perature and low baryon density region in nuclear mat- 
ter phase diagram [l], Q • Several experimental observ- 
ables like the Number-of-Constituent-Quark (NCQ) Scal- 
ing of elliptic flow, jet-quenching, etc. can be explained 
commendably with the appearance of partonic degree 
of freedom in the collisions at RHIC energy regions ^3.]. 
However, due to the complex nature of the relativistic 
nucleus-nucleus reactions, the QGP, if it has been cre- 
ated, escape direct detection. The fact we have not ob- 
served any observation undergoing dramatic change re- 
minds us to perform an energy scan from high energy to 
lower energy to search for the phase boundary and the 
possible end-point of the phase boundary [J. On the 
other hand, in last two decades the heavy ion collisions 
performed at the BEVALAC/LBNL and SIS/GSI 0,0 
were used to produce hot and compressed nuclear mat- 
ter for studying more about the nuclear equation of state 
(EoS) [tI, Si at high baryon density and low tempera- 
ture region of the phase diagram. We have made great 
efforts in studying the nuclear EoS, both theoretically 
and experimentally, but a solid conclusion can hardly be 
made. From this two points for more understanding of 
the nuclear matter phase diagram and EoS at high net- 
baryon density region, it is expected that the Heavy Ion 
Research Facility in Lanzhou china (HIRFL) — Cooler 



Storage Ring (CSR) with a maximal beam kinetic energy 
of 520 MeV/nucleon for heavy nuclei [9|, which focuses 
at the high baryon density region, can make significant 
contribution to those important search. 

The uranium-on-uranium (UU) collision is most suit- 
able for this study. Uranium is the most deformed stable 
nucleus. Representing U as a homogeneous ellipsoid with 
one long (i?;) and two short (i?s) semi-axis, one can re- 
lated their ratio to deformation parameter used in nuclear 
physics 

Ri _ 1 + 4^3 

R, 4-2(5/3^ ^ ' 

For 6u « 0.27, the ratio of the long-axis over short-axis 
is as large as 1.29. For A=238, we will use Ri — 8.4, 
Rs = 6.5 and R = {RlRif/^ = 7.0/m [l^. 

Due to larger A and deformation, the gain in energy 
density for UU over AuAu can reach the factor 1.8 in 
RHIC/BNL [ll|. Because of the deformation, UU colli- 
sions at the same beam energy and impact parameter but 
different orientations are expected to form dense matter 
with different compressions and lifetimes. 

We consider two extreme collisions: the head on colli- 
sions with long-axis on long-axis and short-axis on short- 
axis as tip-tip and body-body collisions, respectively . 
Random collision geometries, which are illustrated in 
Fig.l, lie between them. It is expected that the tip-tip 
collisions can form a higher densities of nuclear matter 
with longer duration than in body-body or the spherical 
nuclei collisions and easy to reach thermal equilibrium at 
the same energy and impact parameter. This is a pow- 
erful tool for studying the physics of large compression, 
high-baryon density and possible phase transition from 
nuclear matter to a new form of matter with partonic 
degree of freedom. 
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FIG. 1: The initial collision geometries of the reaction plane 
in U+U collisions: Z-direction is defined as the beam direc- 
tion. The impact parameter b is along X axis. X-Z plane 
represents the reaction plane. 6p and 9t are the orientation 
of reaction plane for projectile and target, respectively. Thus, 
6p — 6t = O'^ / 90^' corresponded to tip-tip/body-body col- 
lisions respectively. 



For the non-polarized UU collision, target and pro- 
jectile have random orientation at the initial coordi- 
nate space. Several ideas have been developed during 
the past few years to take advantage of the UU colli- 
sions [isl [T^ . uB^ . but no experimental implementation 
has been made. One of the uncertainties is how to select 
collision with the desired geometry such as tip-tip and/or 
body-body. 

In previous Ref [12, [iBl , we have studied systemically 
the stopping power and anisotropy flow in tip-tip and 
body-body UU collisions. In this paper, we apply a rel- 
ativistic transport models, ART [l^, to study the effect 
of different colliding orientations and only focus on the 
central baryon density in order to confirm quantification- 
ally that what initial geometries are the interesting events 
with longtime and high density. This will provide a useful 
help for quantitative analysis in experiment. 

The outline of this paper is as follows. Section II is 
devoted to the centrality dependence of central baryon 
density. In section III, we will discuss the effect with 
different inial collision orientations. In section IV, two 
measurable ways to select out the tip-tip like events will 
been discussed. Finally, a short summary will be given 
in section V. 



II. THE CENTRALITY DEPENDENCE OF 
CENTRAL BARYON DENSITY 

To form a quark-gluon plasma, it is necessary to 
achieve a high local energy or baryon densities in a suf- 
ficiently large volume and for a sufficiently long time so 
that the initial plasma droplets can grow up. The cur- 
rent estimate for the critical baryon and energy densi- 
ties at which the QGP may form are about 5 po and 
2.5 GeV/fm"^, respectively [20|- Under the assumption 



of full nuclear stopping, which is a good approximation 
for heavy nuclei at CSR energy region, as the results of 
energy and baryon number conservation, the maximum 
energy density and maximum baryon density for heavy 
nuclear collisions are 

= 2.5; = 2.5, (2) 

eo Po 

respectively (2l| . 

The full nuclear stopping is reached most possibly in 
central collision. We use the effective centrality which is 
defined with b = b/bmax to replace the impact parameter 
b. Here bmax is the maximum for minimum bias events, 
corresponding to the tip-tip and body-body UU collisions 
are about 13 fm and 17 fm, respectively. 
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FIG. 2: Time evolution of central local baryon density in 
UU tip-tip (solid line) and body-body (dashed line) collisions 
with a beam kinetic energy of 0.52 GeV/u in four different 
centralities. 

Fig. 2 shows the evolution of central baryon density 
in UU tip-tip and body-body collisions at i?boam = 
Q.h2GeV/u at different centralities b using the soft nu- 
clear equation of state with K = 200 MeV. Here, the 
central region is a cube with Ifw? and the lattice esti- 
mate is done for central density [l^. 

As expected, the maximum baryon density decreases 
evidently from central to peripheral UU collisions for 
both tip-tip and body-body. It is remarkable that the 
decline is faster and more obvious in tip-tip than in body- 
body. The orientation of UU collision almost has no effect 
during the early stage of the central collision correspond- 
ing to 6 < 0.3 when the kinetic energy is much higher 
than the potential energy. A maximum baryon density 
of about 3.2 Po a-re reached at about 7 fm/c in both tip- 
tip and body-body collisions at 6 = 0. The matter in 
the high density region (i.e., with p/po ^ 2.5 that may 
occur the full nuclear stopping) lasts for about 13 fm/c 
and 25 fm/c at body-body and tip-tip coUisions, respec- 
tively. This means that the high density phase in the 
tip-tip collisions lasts about two times longer than that 
in the body-body collisions. But this discrepancy of high 
density phase lifetime fades away rapidly from 6 = to 
b — 0.3, even the high density phase vanishes for tip-tip 
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collisions at about b = 0.6. The higher compression and 
longer passage time render the central tip-tip UU colli- 
sions as the most probable candidates to form the QGP 
and to study nuclear EoS in high density condition. 

One should take note of that the tip-tip and body- 
body UU collisions can reached a near maximum baryon 
density but exits a different lifetime of high density mat- 
ter for central events (See Fig2(a)). In order to confirm 
what events are the interesting events with extended life- 
time, this discrepancy of tip-tip and body-body must be 
displayed quantificationally. We define a new variable — 
the baryon aggradation strength as followed: 



5* 



t' 



pdt, 



(3) 



Here, t and t are the threshold and end time of the 
central high density matter. The centrality dependence 
of baryon aggradation strength is showed in Fig. 3. One 
can see clearly that the strength of tip-tip is higher than 
body-body at & < 0.4. The curve for tip-tip is flat from 
6 = to 6 = 0.2 but has a sharp drop a,t b — 0.2. Oth- 
erwise, the high density matter vanishes for tip-tip but 
still exist for body-body at b > 0.7. Therefore, as a con- 
servative estimate, the tip-tip events with b < 0.2, corre- 
sponding to 6 < 2.6 fm and the value of S" > 70, are the 
interesting events with extended high density lifetime. 

In Fig. 3, we also indicate the result of using the forward 
(i.e. 9 < 15°) neutron multiplicity 7V„ to replace the 
normalized impact parameter b because the multiplicity 
of forward neutron which can been measured easily in 
experiment at CSR and have a linear relation with the 
impact parameter One can see clearly that the 

multiplicity of forward neutron in tip-tip is smaller than 
in body-body at 6 < 0.2. Hence, we can separate easily 
the tip-tip and body-body events only by the cut of Nn < 
30 when most tip-tip central events will survive. 



III. THE ORIENTATION DEPENDENCE OF 
CENTRAL BARYON DENSITY 

From Fig.l., one can see clearly that it depends 
strongly on five parameters to confirm the initial geome- 
tries of UU collision. They are the impact parameter b, 
the initial orientation of projectile and target in reaction 
plane 9p and 9t, the inial orientation of projectile and 
target in transverse plane (/)p and (/)t, respectively. 

In section II, we gained a simple estimate that the 
interesting events with extended high density lifetime re- 
late to 6 < 2.6/to. Based on this point, we will discuss 
the orientation dependence of central baryon density un- 
der b < 2.6/to in this section. 

Fig. 4 (a) and (b) are the orientation of initial trans- 
verse plane dependence on baryon aggradation strength 
at 6 = fm and b = 2.6 fm, respectively. Here, we only 
need to take into account the relative orientation in ini- 
tial transverse plane because of the symmetry of initial 
overlap region (interaction region) in shape. 
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FIG. 4: The relative orientation of projectile and target in 
the initial transverse plane dependence on baryon aggradation 
strength at 6 = fm (a)and 6 — 2.6 fm (b), respectively. 
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FIG. 3: Left: Baryon aggradation strength versus normal- 
ized impact parameter for UU tip-tip (circular hollow dot) 
and body-body (triangular solid dot) collisions. Right: The 
forward neutron multiplicity versus normalized impact pa- 
rameter for UU tip-tip (dashed line) and body-body (real 
line)collisions at a beam kinetic energy of 0.52GeV/u. 



From Fig. 4, we can see that the baryon aggradation 
strength is nearly invariable when the orientations of pro- 
jectile and target (0p, 9t) in the reaction plane are fixed 
but the orientations of projectile and target in the trans- 
verse plane {(j)p, (jir) are changed from 0° to 90° at both 
6 = and b = 2.6/to. Here, we only consider the ex- 
treme condition that the orientations of projectile and 
target are equal, i.e. 9p = 9^- There is an express drop 
is showed from O" to 30° for 9p or 9t when the rela- 
tive orientation in transverse plane Ac/) is fixed. This 
phenomenon illuminates both 9p and 9t should are less 
than 30° for the tip-tip like event. 

Due to inessential (pp and (/)t, fig.5(a) are the initial 
orientation dependence on baryon aggradation strength 
at = fm with (f)p — (j>T = 0°. Hence, we only need 
to consider two extreme geometries for full overlap in 
transverse plane. Others are between them. The curve 
indicates the discrepancy is tiny for 9p = ±9^ at b = 
0/to. It is remarkable that a sharp drop at 9p = \9t\ — 
20° is seen. If one requires S > 70, the orientations in the 
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those events with biggish impact parameter b. 
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FIG. 5: The initial orientation of projectile and target in 
reaction plane dependence on baryon aggradation strength in 
UU collisions at Eb = 0.52GeV/u. Left window (a) with same 
initial orientation 9p — \9t\ &t b — Ofm. Right window (b) 
with fixed Or = 20° at b=0 and b = 2.6/m. 



FIG. 6: Left window (a) is the forward neutron multiplicity 
A^„ distribution with different selective impact parameters b. 
Right window (b) is the impact parameter b distribution with 
different selective forward neutron multiplicity in random UU 
collisions at Et = 0.52GeV/u. 



initial transverse plane should be less than 20° for both 
projectile and target. Fig. 5 (b) shows 9p dependence on 
baryon aggradation strength with fixed 6t = 20" at both 
b = Ofm and b = 2.6fm. It emphasizes again that 9p is 
less than 20" if 5* > 70. 

In a short word, those events with extended high den- 
sity phase in UU collisions are the tip-tip like events 
which are with the initial orientations dp < 20° and 
6t < 20° in the transverse plane and b < 2.6/m. 



IV. TWO MEASURES FOR SELECTING OUT 
THE TIP-TIP LIKE EVENTS IN EXPERIMENT 

Due to the effect of initial orientation of projectile and 
target for deformed nucleus, one key point to program 
UU collisions in experiment is how to select out those 
collisions with the desired geometry such as tip-tip like 
events. 

In view of the discrepancy of forward neutron multi- 
plicity Nn dependence on impact parameter b between 
tip-tip and body-body UU collision showed in fig. 3, a 
available try is associated spontaneously. Fig. 6(a) shows 
the forward neutron multiplicity iV„ distribution with dif- 
ferent selective impact parameters b in random UU colli- 
sions. Here, About 6.8 millions events with b < 6fm are 
created for our simulation. Obviously, it can effectively 
reject those events with b > Afm by selecting Nn < 50. 
Fig. 6(b) shows impact parameter b distribution with dif- 
ferent selective forward neutron multiplicity. If one se- 
lects the cut with Nn < 40, the percentage of those events 
with b < 2.6/m will been enhanced. It is worth to men- 
tion that the percentage can been farther enhanced with 
the cut Nn < 30, but large numbers of tip-tip like events 
with b < 2.6/to also will been rejected. Thereby, Nn is a 
good measure in experiment as a trigger for fast rejecting 



As a incidental available result, fig. 7 shows the initial 
orientation {0p and 9t) distribution of projectile and tar- 
get in reaction plane with the cut iV„ < 40 and Nn < 30 
for random UU collisions. One can see the most events 
is tip-tip hke events with both 9p < 20° and 6t < 20°. 
It is satisfying by the reducing Nn to gather more pure 
tip-tip like events. 



(a) N <40 



(b) N <30 





FIG. 7: The initial orientation {6p and 6t) of projectile and 
target in reaction plane when cut A'„ < 40 (a) and Nn < 30 
(b) in random UU collisions at E}, = 0.52GeV/u. 

In addition, it is well-know that an obvious long dura- 
tion will load possible thermal equilibrium of high density 
matter which has been obtained in tip-tip central colli- 
sions. The degree of thermalization can be measured by 
the ratio R of transverse to longitudinal momenta in low 
and intermediate energies heavy ion collisions [l7j . Its 
expression is 



2E. \P^■. 



(4) 



Here, i is the serial number of particles, the sum involves 
all particles, \pit\ and \piz\ are the total absolute value 
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of nucleon transverse and longitudinal momentum in the 
c.m.s, respectively. It's a multi-particle observable on an 
event-by-event basis, which for an isotropic distribution 
is unity. Then i? = 1 is a necessary, although not is a 
sufficient, condition for thermal equilibrium. It is also a 
measure of the stopping power [18| . 



thrown away. 

Fig. 9 shows the initial orientation (dp and 9t) distri- 



1.2- 




10' 



10^ 



10' 



Random 


1 : 

(b) 




L b<2.6fm - 


I|8pep|>70° 
















^ / \ 





2 4 6 
Impact Parameter b (fm) 



0.6 0.8 1 1.2 1.4 
Stopping Power R 



FIG. 8: Left: the nuclear stopping power R as a function of 
the impact parameter b in different orientation UU collisions 
at iSi, = 0.52GeV/u. Right: the nuclear stopping power R 
distribution with b < 2.6 fm in different orientation UU colli- 
sions at Ei = 0.52GeV/u. 

Fig. 8 (a) shows that the impact parameter b depen- 
dence on nuclear stopping power R in different orien- 
tation UU collisions. Here, Only charge particles are 
considered. For three samples with different orientation 
collisions, a nice linear relation can been seen clearly be- 
tween stopping power R and impact parameter b. For 
the most central events, tip-tip like collisions {9p^T < 20) 
have a maximal value for R among three samples, i.e. 
maximal most strong stopping power. When b increase, 
R decrease. Otherwise, that stopping power R in tip- 
tip like collisions is large than in body-body like events 
{{Op,T > 70) at 6 < 4/m. Even, R reaches to ~ 1.2 in tip- 
tip like events but only ^ 0.92 in body-body like events 
and ~ 1.02 in random collisions at 6 = 0/m. Thus, we 
can gather the center events by selecting R, for instance 
the cut with R > I corresponding to & ~ Ifm for random 
and b ~ 2/m for tip-tip like UU collisions can get rid of 
body-body like events. 

This remarkable discrepancy prompts us that the cen- 
tral tip-tip like events maybe have a distinct distribution 
of stopping power R. Fig. 8(b) shows stopping power R 
distribution with b < 2.6 fm for three samples with dif- 
ferent orientation UU collisions. It can be see clearly 
that the body-body like events can been rejected com- 
mendably from random events at 6 < 2.6 fm by selecting 
R > 1. If one expects higher percentage of tip-tip like 
events in the sub-sample, i? > 1.1 can been selected, but 
at the same time vast tip-tip like events also will been 



(a) R>1 



(b) R>1.1 




FIG. 9: The initial orientation {6p and 9t) of projectile and 
target in reaction plane when cut i? > 1 (a) and i? > 1.1 (b) 
in random UU collisions at Et = 0.52GeV/u. 



bution of projectile and target in reaction plane when 
cut i? > 1 (a) and R > 1.1 (b) in random UU collisions. 
It can been seen very clearly that most events have a 
small orientation, which events are closer to tip-tip like 
events, by selecting the nuclear stopping power R > 1 
and R > 1.1. When R > 1.1, the events with larger 
orientation are rejected ulteriorly. In other words, it can 
enhance the purity of tip-tip like events and minify the 
background when R increases. Therefore, with the nu- 
clear stopping power R, it is possible to select out the 
tip-tip like collisions in experiment. This will allow us 
to study the nuclear collisions at higher density with a 
considerable long duration. 



V. SUMMARY 

Performing UU collisions are a good tool to search for 
QGP, the phase boundary and the possible end-point of 
the phase boundary. Those most interesting events (i.e. 
tip-tip like) is with small initial orientations (< 20*^) in re- 
action plane and small impact parameter (< 2.6 fm). For 
the non-polarized UU collisions, we have developed two 
measures allow us to select out those tip-tip like events in 
experiment. The forward neutron multiplicity as a trig- 
ger and nuclear stopping power is suited for fast on-line 
and off-line physical analysis, respectively. 
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